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The hard photon productions from bremsstrahlung and annihilation with scattering that arise at 
two loop level are estimated from a chemically non-equilibrated quark gluon plasma using the frame 
work of Hard Thermal Loop (HTL) resummed effective thermal field theory. The rate of photon 
production is suppressed due to unsaturated phase space compared to it's equilibrium counterpart. 
However, the suppression is relatively weaker than expected from the kinetic theory due to an addi- 
tional collinear enhancement arising from the decrease in thermal quark mass. For an unsaturated 
plasma, unlike the effective one loop case, the reduction in the effective two loop processes is found 
to be independent of gluon fugacity, but strongly depends on quark and anti-quark fugacities. It 
is also found that, since the phase space suppression is highest for annihilation with scattering, the 
photon production is entirely dominated by bremsstrahlung mechanism at all energies. This is to 
be contrasted with the case of the equilibrated plasma where annihilation with scattering dominates 
the photon production particularly at higher energies. 



PACS number(s): fl.IO.Wx, I2.38.Mh 
I. INTRODUCTION 

The study of single photon production at relativistic 
heavy ion collisions has gained momentum in recent years 
due to availability of experimental data from CERN, SPS 
and also the data expected shortly from the RHIC exper- 
iments at BNL |I] . Assuming the formation of a quark 
gluon plasma (QGP), the theoretical studies utilize the 
effective field theoretical formulation with hard thermal 
loop (HTL) resummation technique |f^-[Tl[| to calculate 
the imaginary part of the photon self energy and hence 
the photon rate. An important aspect in this approach 
is to distinguish hard momentum of order T from soft 
momentum of order gT where g is the QCD coupling 
constant. The propagation of soft momentum is con- 
nected with infrared divergences in loops and therefore 
the propagators need to be dressed to get finite result. 
According to Braaten and Pisarski j^J^] , for soft momen- 
tum, instead of using bare propagators and vertices, ef- 
fective propagators and vertices should be used. This 
method has been adopted for the calculation of the rate 
of soft dilepton production [Of and hard real photon pro- 
duction due to annihilation [qq — > 7<?) and QCD Comp- 
ton (qg — > 57, qg — > qj) processes from a quark matter 
at one loop level in the effective theory Jl^Jl^]. In this 
method, a cutoff parameter (gT < k c <T) is introduced 
to distinguish between the soft and the hard quark mo- 
mentum circulating in the loop. For hard real photon 
production, it is sufficient to use an effective propaga- 
tor (summed over successive one loop insertions) for one 
of the quark loops carrying momentum below the cutoff 
while the other loops and vertices can remain undressed. 
Above the cutoff, bare propagators and vertices can be 
used and a loop correction must be inserted on the hard 
propagator. When adding the soft and the hard contri- 
butions, the cutoff dependence cancels out. 

Alternatively, the rate of photon production can also 



be estimated on the basis of relativistic kinetic theory 
where the integration is carried out over a phase space 
volume multiplied by the square of the reaction ampli- 
tude and the appropriate distribution functions for ini- 
tial and final states. In this approach also a cutoff k c 
for integration can be introduced so that the soft part 
that involves divergence can be treated separately. It is 
interesting to note that the total photon production rate 
can also be estimated directly from the hard part using 
a lower cutoff parameter for integration equal to twice 
the thermal quark mass (k% — 2m^) |l3). This approach 
has been extended to calculate photon production from 
non-equilibrium plasma |L^JL6[ ] with the use of thermal 
quark mass and parton distribution functions appropri- 
ate for a non-equilibrium situation |L7| . In case of a non- 
equilibrated plasma, additional contribution is expected 
from the pinch singularity |n] Q9L However, it is shown 
in ]17t that pinch contribution at soft momentum scale 
is subleading with respect to the dominant HTL contri- 
bution. Similarly, for the hard scale, pinch singularity is 
absent due to restricted kinematics. 

It may be mentioned here that in the above cutoff 
method, a bare gluon propagator has been used even if 
the cutoff does not constrain the gluon to be hard. Allow- 
ing the gluon to be soft, leads to new physical processes 
that may contribute to the hard photo n p roduction. Re- 
cently, it is shown by Aurenche et al. |2(J,|2l| that signif- 
icant contribution comes from the bremsstrahlung and 
a new process called annihilation with scattering (AWS) 
that arise at two loop level in the effective theory due 
to space like soft gluon exchange. Their results can also 
be written in a way that separates the phase space from 
the amplitude of the process producing the photon. The 
magnitude of the amplitude usually becomes less when 
the number of loops increases. However, it is found that, 
within the effective theory, the phase space contribution 
at one loop level turns out to be smaller than the two 
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loop due to kinematical constraints. Both effects com- 
pensate so that two loop diagrams in the effective theory 
also contribute at the dominant level. 

The work of Aurenche et al. (2(],|2lJ assumes the plasma 
to be in equilibrium at temperature T. However, the 
rate of photon production may be affected significantly if 
the phase space remains unsaturated. The present work 
extends the formulation of Aurenche et al. to the non- 
equilibrium QGP. We estimate the photon production 
from bremsstrahlung and AWS processes for a chemically 
unsaturated quark gluon plasma. We restrict to the re- 
gion of Landau damping part (L 2 < where L is the 
gluon four momentum), whereas the region I? > for 
hard gluon exchange has been included in the one loop 
calculations in the effective theory. In a subsequent work 
p2| , p3| , Aurenche et al. have shown that even the higher 
order contribution can not be ignored for real photon pro- 
duction indicating that the thermal real photon produc- 
tion in QGP is a non-perturbative mechanism. On the 
other hand, in such situation the applicability of HTL 
resummation technique which is based on perturbative 
approaches, may become questionable However, the 
purpose of the present work is not to go into the above as- 
pect in detail. Here, we only focus on the bremsstrahlung 
and the AWS photon production from a chemically un- 
saturated quark gluon plasma which is of significance at 
RHIC and LHC energies |§-|||. 

First we consider the bremsstrahlung and AWS photon 
productions from an equilibrated QGP p4| . We draw 
similar conclusions as that of previous work that, 
within the framework of effective theory, the two loop 
contributions particularly due to AWS, compete with one 
loop contributions at all energies. For a chemically non- 
equilibrated plasma, since the phase space is unsaturated, 
the photon productions both at effective one and two 
loops level are suppressed as compared to the equilibrated 
case. Since the thermal quark mass decreases with fugac- 
ities, the collinear enhancement for the effective two loop 
processes also goes up. As a consequence, the suppres- 
sion at the effective two loop level has been found to be 
weaker than expected. Interestingly, the above suppres- 
sion is independent of the gluon fugacity and depends 
only on unsaturated quark and anti-quark distribution 
functions. Further, it is noticed that the suppression for 
the AWS process is the highest and the photon produc- 
tion is entirely dominated by the bremsstrahlung mech- 
anism particularly when the plasma is strongly unsatu- 
rated. This is contrary to the case of equilibrium sit- 
uation where AWS is the dominant mechanism of pho- 
ton production at higher energies. This result can be 
qualitatively understood on the basis of relativistic ki- 
netic theory where the AWS process has more incoming 
fermion lines compared to the bremsstrahlung process. 
We may mention here that based on kinetic theory ar- 
gument and using the equilibrium results, an extension 
to non-equilibrium QGP has been discussed by Mustafa 
et al. p5[ . However, their results are not in agreement 
with the present findings which are derived using the for- 



malisms given in | ]20L|21[ . 

The paper is organized as follows. We begin with the 
description of an unsaturated plasma with a brief review 
of the photon production both at one and two loop level 
in the effective theory in section II. In section III, we 
evaluate the bremsstrahlung and the AWS photon pro- 
duction from an unsaturated quark gluon plasma. We 
show that the imaginary part of the self energy can be 
written in a form which separates the amplitude of the 
reaction from the phase space so that the use of kinetic 
theory can be justified for non-equilibrium plasma. We 
discuss the results of photon production both at one and 
two loop levels in the effective theory in section IV. Fi- 
nally, the conclusion and summary are presented in sec- 
tion V. For effective two loop calculation, we follow the 
Retarded/ Advanced (RA) formalism [ ^6p^ ] where the 
propagators remain same as zero temperature field theory 
while the vertices are redefined which include nonequi- 
librium distribution functions. It is shown in the ap- 
pendix that the redefined vertex has the same form as 
that of equilibrium case except that the distribution func- 
tions need to be defined appropriately to represent a non- 
equilibrium phenomena. 



II. GENERAL FORMALISM 

We consider a thermalizcd plasma of quarks and glu- 
ons expected to be formed during the collisions of two 
heavy ions at relativistic energies. However, at RHIC and 
LHC energies, several perturbative-inspired QCD models 
p8| , [39f predict the formation of an unsaturated plasma 
with high gluon content j p9j . Such a plasma will attain 
thermal equilibrium in a short time to w 0.3 — 0.7 fm, but 
will remain far from chemical equilibrium [ p5[ . Since the 
initial plasma is gluon dominated, more quark and anti- 
quark pairs will be needed in order to achieve chemical 
equilibration. The dynamical evolution of the plasma un- 
dergoing chemical equilibration was studied initially by 
Biro et al. [ p5[ and subsequently by many others [jl6 33 
by solving the hydrodynamical equations along with a set 
of rate equations governing chemical equilibrations. In 
this work, we do not consider the hydrodynamical evolu- 
tion of the plasma, rather calculate only the static rate 
of real photon production from an unsaturated plasma. 
Further, we also assume an ideal situation where the 
plasma is baryon free. 

A chemically non-equilibrated but thermally equili- 
brated plasma can be described by the Juttner distribu- 
tion function for quarks (anti-quarks) and gluons, given 
by, 
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The fugacity factor A;, (i = q,q,g) is related to the 
chemical potential fii as Xi — e p " i l T . The plus and minus 
signs are meant for the Fermi-Dirac and Bose-Einstein 
distributions respectively. At equilibrium, chemical po- 
tential vanishes and Xi — > 1. The distribution func- 
tions can also be factorized in an approximate way as 
rii = Xin eq where n eq is the distribution function at equi- 
librium. In this representation, Xi gives a measure of 
deviation from the corresponding Fermi-Dirac or Bose- 
Einstein distributions. In case of a baryon rich plasma, 
the quark and the anti-quark distribution functions can 
also be described by the same Juttner distribution func- 
tions. However, the chemical potential will have an addi- 
tional component to account for the finite baryon density 
of the plasma (See Refs. [[H||3| and the appendix A for 
detail ). In the present case, since the plasma is baryon 
free, the fugacity and the chemical potential both refer 
only to the unsaturation properties of the plasma. 

The thermal photon production rate from such a 
plasma can be related to the retarded polarization ten- 
sor of the photon using thermal field theory [^(J. For 
real photons, this relation gives the number of photons 
emitted per unit time and per unit volume of the plasma 
as 
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where E = qo is the energy of the emitted photon and 
II is the retarded self-energy at finite T. This relation is 
valid only at first order in the QED coupling constant a 
but is true for all orders of the strong coupling constant 
a s , as it has been assumed that the produced photons 
emerged from the matter without further scattering. It 
may be mentioned here that the emitted photon of en- 
ergy E is not the part of the heat bath and therefore, has 
no distribution of any form. Thus n B (E) in Eq. (|4|) is 
just a prefactor which happens to coincide with the form 
of Bose-Einstein distribution only when the plasma is in 
full equilibrium. The photon production due to QCD 
Compton and annihilation processes can be estimated 
by evaluating the photon self energy from the diagrams 
shown in Figure 1. The imaginary part of the self en- 
ergy can be obtained by cutting the above diagram using 
thermal cutting rule It may be mentioned here 

that the rate of photon production evaluated from imag- 
inary part of photon self energy with some finite order 
of loop expansion is equivalent to the relativistic kinetic 
theory [jl3| . The self energy calculated upto N loop level 
for m particles — > n particles + 7, is equivalent to the ki- 
netic theory estimates from all reactions consistent with 
m + n < N+ 1. 





FIG. 1. Diagrams contributing to the photon self energy for 
the Compton and annihilation processes when intermediate 
quark momentum is (a) hard (b) soft 

Similarly, the photon production from bremsstrahlung 
and annihilation with scattering (AWS) processes can be 
obtained from the two loop effective diagrams as shown 
in Figure 2 (notations are same as given in |20| , |2l|| ) . 




FIG. 2. (a) Vertex and (b) Self diagram for two loop self 
energy in the effective theory 

In section III, we estimate the above self energy (at ef- 
fective two loop level) more explicitly for an unsaturated 
plasma as defined above. 



III. PHOTON PRODUCTION AT TWO LOOP 
LEVEL 

In the two loop level of effective theory the Feynmann 
diagrams need to be considered are shown in Figure 
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g. The imaginary part of the photon self energy can 
be expressed as a sum over the possible cuts through 
the effective two loop diagrams. The physical processes 
bremsstrahlung and quark anti-quark annihilation are 
obtained by cutting through the effective gluon propa- 
gator. To obtain the contribution of bremsstrahlung and 
AWS, the same order of magnitude as of contribution 
from Compton and annihilation processes, the quark mo- 
mentum circulating in the loop should be hard. As a re- 
sult, all the vertices and propagators can be used as bare 
one except for the gluon propagator since the gluon can 
be soft. Recall that only Landau Damping part (L 2 < 0) 
gives bremsstrahlung and AWS, whereas the L 2 > part 
gives Compton and annihilation process which has been 
already included in the effective one loop level calcula- 
tions for hard gluon exchange. However, due to phase 
space restriction, the Landau damping is the dominant 
mechanism if L is soft pTt . 

Figure ^| shows the relevant cuts and circling required 
to evaluate the self energy using thermal cutting valid 
for RA formalism |3(||37],[l3]] . In the following, we extend 
the formalism of Aurenche et al. [^,|T) to the chemi- 
cally non-equilibrated situation. An important aspect in 
the RA formalism is the redefined vertices which con- 
tain the distribution functions. We have derived them 
in appendix A for a more general situation where the 
plasma is chemically unsaturated as well as has non zero 
baryo-chemical potential. However, in the present work, 
we consider only a baryon free chemically unsaturated 
plasma. The vertex functions corresponding to Figure 
2(a) are given by 

e AAR (R, -P, -Q) = e[n F (r„) - n F ( Po )} 
g AAR (R + L, —L, -R) = g[n B (l ) + n F (r + I )] 
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-g ARR (-P - L,P,L) 



= -9 (5) 

The above vertex functions are in the same form as the 
equilibrated case except the distribution functions should 
contain appropriate chemical potential. In the above ex- 
ample, the chemical potentials associated with R and P 
lines are /j, q where as it vanishes for Q and L lines. There- 
fore, using the above definitions, the vertex (Figure §(a)) 
and self (Figure |^(b)) diagrams of the imaginary part of 
the photon self energy can be expressed as 
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According to the cutting rule valid for RA formalism jl3j , 
the fermion cut propagators are given by 

„RA , v A , 

S (P) = S (P) 

AR , A , , 

S-(P) = -S (P) 

s a -(p) = s r (p)-s a (p) 

S RA (P)^S A (P)-S R (P). (8) 
We denote the fermion propagator: 



S ' (P) = IfS • (P) with P ee ( Po , ^/p 2 + Mlp) (9) 



(10) 



P~ ± i Po e P 2 ~ M lo ± Woe ' 
Similarly, boson cut propagator 

D A f(L)^D R a (L)-D A a (L) . (11) 

The effective gluon propagator in a linear covariant gauge 
is given by 

~~R,A , . T , . . R.A , . L , , . R,A , 

-D pa (L) = P pa (L)A T - (L) + P pa (L)A L (L) 
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Here P pa are the usual transverse and longitudinal pro- 
jectors in linear covariant gauges and 

c 2 x(l - x 2 ) , ( x + 1 



n T (L) ee 3m z 



U L (L)EE3m 2 (l-x 2 ) 



■In 



1 



X 


^ + 1V 








K x-l)\ 



(13) 
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where, x — Iq/1. 

Using the expression of the vertices in Eq. (||) and the 
cut propagators, the vertex diagram of photon self energy 
can be simplified as 
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Similarly, the photon self energy for the self diagram can 
be expressed as in Ref. |2l| with appropriate distribution 
functions 
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where e is the electric charge of the quark which depends 
on its flavor. The trace in Eq. ( p"5| ) and Eq. (plf) are given 
by, 



Traced Uf = Tr[ 7 ^ 7M (#) 7 "(# + R] . 



(17) 



The factor S(R), S(P + L) without any R and A super- 
script denotes the principal part of the propagator i.e., 



S(R) = 
S(P + L) = 
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The difference between the retarded and the advanced 
gluon propagators, in both the vertex (Eq. |l5|) and self 
(Eq. |l^) part of the photon self-energy is known as spec- 
tral function given by 
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The properties of these spectral functions depend upon 
the analytic structure of the gluon propagator. For L 2 < 
(i.e. \x\ < 1) region, in which we are interested about, 
the self-energies ±1 T L acquire an imaginary part due to 
the logarithm and the corresponding expression of the 
spectral functions are given by, 
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Here the imaginary part of retarded self energies are 
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Adding the contributions from these two diagrams and 
plugging into the expression (Eq.(|4|)), the photon produc- 
tion from the effective two loop level can be evaluated. It 
should be noticed that the expression for the self energy 
is same as used in |21 , except for the non-equilibrium dis- 
tribution functions n F and thermal masses which contain 
chemical potentials. Like in one loop case, here also we 



ignore the pinch contribution and assume that the effec- 
tive quark and gluon propagators are still given by their 
equilibrium counter part with use of asymptotic thermal 
quark mass = 2m 2 q defined by 
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3tt 2 



p dp [2n g + (n q + rig)} 



(22) 



which can be shown to be (A 9 + -j-)g 2 T 2 /9 for factorized 
distributions. The thermal gluon mass m 2 appropriate 
for non-equilibrium plasma Q,^6| is given by 



^2 / PMGn g 



Nf(n q +n q )} 



(23) 



For a factorised distribution the m 2 can be written as 

(3A 9 + ^j-X q )g 2 T 2 /9. The kinematic conditions restrict 
the phase space for the physical process into three regions 
as shown in Figure [§[ 




(b) (c) 

FIG. 3. Physical processes appear by cutting the two loop 
diagram of Figure 2 for L 2 < (a) Region I: (b) Region II (c) 
Region III 

Region I p < and r + l Q < 
Region II po < and ro + Iq > 
Region III p > and r + l > 

Region I and III corresponds to bremsstrahlung from 
anti-quark and quark respectively. The region II cor- 
responds to AWS processes. The physical processes with 
the lower line in the Fig. || replaced with quark, anti- 
quark have also been included. It may be mentioned here 
that region I and region III will give same contribution as 
long as quark and anti-quark distribution functions are 
same (i.e. for a baryon free plasma). We will study the 
contribution from region III and multiply it by a factor 2 
to get the total bremsstrahlung photon yield. Similarly, 
we will discuss about the region II for AWS process. 



A. Bremsstrahlung 

The two points which need to be addressed here are the 
collinear limit and the enhancement due to collinear sin- 
gularity. It is important to notice that due to the factors 
[(R 2 - Mi]((P ' + L) 2 - Af^)]" 1 which appear in the inte- 
gral Eq. (ffq) is responsible for the collinear divergence. 
Although we will carry out the complete integration, even 
at the qualitative level, we can understand the nature of 
the divergence from the integral of the type 
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du 



2tt 



(24) 



where (u = 1 — cos9) is the angle between p and cf and 
<p is the azimuthal angle between q and / when projected 
on a plane orthogonal to r. For hard photon production, 
the following approximations can be used 
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In the above, (b — a) ~ L 2 jp 2 gives a measure of the 
distance between the two poles. Since the two poles are 
very near by (for soft gluon exchange), the above integral 
is of the order 
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p(p 2 q 2 )a 
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for a > L 2 /p 2 or Ml > L 2 . Note that this enhancement 
factor which is associated with the smallness of the angle 
of emission is same irrespective of whether the emitted 
photon is soft or hard. As a consequence, the integral is 
enhanced by a factor of order p 2 /M^ ~ l/(g 2 X g ) if the 
plasma is strongly gluon dominated i.e. X q < X g <C 1. 
The above enhancement is larger by a factor of 1 as 
compared to the equilibrated case. Although the above 
argument is qualitative, it remains valid even after all 
the integrations are carried out (see Figure 4 in section 
IV). Now, proceeding as before |2l| under collinear ap- 
proximation, the expression for the imaginary part of self 
energy can be written as 

Im It 4 V 

~ (-^r^eV^ [n F (p)-n F {p + E)\ 

x(p 2 + (p + Ef) dp I l A dl ( dx [ du n B {lx) 
Jo J-i Jo 

xp T , L (Z, lx)(l - x 2 f(l - u')- 1 ' 2 {AM 2 00 + l 2 (l - xVr 1 

(27) 

where u' = ~8r 2 u/L 2 and the symbol (— 1)t denotes an 
extra minus sign in the transverse contribution. Here we 
have introduced some cut-offs p* and I* at a scale inter- 
mediate between gT and T, where we assume r to be hard 
and I to be negligible compared to T. We have ignored 
the factor n F (ro + Iq) since it is much smaller compared 
to the Bose distribution n B (lo) (where Iq — Ix) particu- 
larly for unsaturated QGP. As a result, the integral over 
p becomes independent of integrals over L. The Bose- 
distribution is also approximated to n B (lx) ~ j-. This 
assumption considerably simplifies the numerical evalu- 
ation of the self energy which can be expressed in terms 
of the dimensionless constants Jt and Jl, 
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where 
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Note that we do not assign any chemical potential to the 
gluon line L leading to the approximation £• rather than 
\ g -f£ (see subsection D for more discussion). The func- 
tions Jt, Jl depend on the thermal mass ratio m 2 /Ml 
and I* /Moo. Notice that the above expressions for Jt 
and Jl are same as in Ref. except that the ther- 
mal masses now depend on the chemical potentials or 
fugacities. It has been shown in (2Tjj2l[] that for equili- 
brated plasma, taking w* — > oo introduces a negligible 
contribution to the integration of Jt,l- This argument 
is also valid for chemically non-equilibrated plasma since 
the value of w* increases with decreasing fugacity. Thus 
extrapolating the upper limit of the integration w* to oo 
introduces smaller error as compared to the equilibrated 
plasma and can be neglected. Finally, the imaginary part 
of the self energy can be written as 

^NC F a 2 T 
~ -^4~ e 9 (Jt - Jl)^ 

x [n F (p)-n F (p + E)](p 2 + (p + E) 2 ) dp (29) 
Jo 

The total photon production rate for bremsstrahlung pro- 
cess can be evaluated as 
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dp (p 2 + (p + Ef)[n w (p) -n F (p + E)] (30) 



where e/ is the electric charge of the quark flavor / in 
units of electron charge and n B (E) — [exp(_E/T) — 
is the prefactor independent of any chemical potential 
for the bremsstrahlung process. The Jt and Jl inte- 
grals depend only on the thermal mass ratio m^/Ml 
and are insensitive to the chemical potential or fugac- 
ity. Since the prefactor ub(E) is independent of fi for 
bremsstrahlung, the chemical unsaturation effect enters 
only through the quark distribution functions uf{p) in 
the p integral. Therefore, for factorized distribution 
functions, the bremsstrahlung contribution from non- 
equilibrated plasma is suppressed by a factor of \ q as 
compared to the equilibrium case. 



G 



B. Annihilation with scattering 

In case of qq annihilation with scattering (AWS), one 
should consider the region II where po < 0. The pro- 
cedure to estimate AWS rate is nearly identical to that 
of bremsstrahlung except for the exchanges u — > v and 
p — > — p where v = 1+cos 9. In this case also the enhance- 
ment mechanism remains same as before. Finally, with 
the replacement of p by — p, the p integral in Eq. ( |27| ) 
can be written as 



2E 



dR 
d?q~ 



NCpOiOts 



T 



AWS 



xn B (E) / dp(p 2 + (E~p) 2 )[n F (-p)-n F (E~p)] (31) 
Jo 

where n F (—p) = 1 — n F {p). For equilibrated plasma, 
the contribution from [n F {— p) — n F {E — p)\ is assumed 
~ 1. This approximation is also valid in case of non- 
equilibrated plasma since the distribution functions are 
smaller by a factor of \ q . The unsaturation effect appears 
through the prefactor n B (E) which goes as [exp(£'/T) — 
1] _1 A^ under factorized approximation (see the discus- 
sions in appendix A). Therefore, the AWS photon pro- 
duction is suppressed by a factor of X 2 as compared to 
its equilibrium counterpart. 

Finally, we would like to end this sub-section with the 
remark that both bremsstrahlung and AWS photon pro- 
ductions for non-equilibrated plasma are suppressed by a 
factor of \ q and X 2 respectively due to unsaturated quark 
and anti-quark distribution functions. The suppression 
due to unsaturated gluon distribution function seems to 
get compensated by the additional enhancement caused 
by collinear singularity. 



C. Factorisation and Kinetic theory 

In this section, we establish a connection between 
the field theoretical formalism as given in previous sub- 
section and relativistic kinetic theory which can also be 
used to estimate the photon production rate under semi- 
classical approximation . Consider the case of pho- 
ton production due to bremsstrahlung. In the kinetic the- 
ory approach, the photon production rate can be evalu- 
ated by integrating the amplitude squared of the process 
over the phase space of unobserved particles and given 
by 




x2nS(P 2 - M 2 ^ 2nS((R + Lf~M 2 00 ) 

x2nS(K 2 - Ml) 2ttS((K + L) 2 - A&) 

xn F (r + l ) n F (k ) [l-n F (p )] [1 - n F (k + Z )] 

(32) 

Here we have only considered the amplitude for the 
bremsstrahlung process where quark has been scattered 
from another quark. In order to get total photon produc- 
tion rate from bremsstrahlung, the processes involving 
quark scattered from a gluon or an anti-quark also have 
to be considered. 

In order to establish a connection to field theory, we 
now look for various statistical factors that appear in 
Eq.(^||). Let us consider the product of all the distri- 
bution functions that appear in the calculation of self 
energy [see Eq.(|) and Eq.(|^)] 

n B (qo)[n B (Z ) + n F (r + l )] 

x [n p (r ) - n F (p )] [n p (k + Z ) - n F (k )] 

= n F (r a + l a )[l - n F (p )]n F (k )[l - n F (k + l )] (33) 

The pre-factor n B (qo) (go = E) comes from Eq.(^) while 
the second and third factors come from the vertex func- 



tions g 



a An 



and 



^AAR 



Although not explicit in Eq.([15| 



the last factor appears due to the hard thermal quark 
loop contribution to the gluon self energy 01 . In the last 
factor, n F will be replaced by n g when the quark scat- 
tering from gluon is considered. For the above identity 
to be valid, the distribution functions should be defined 
properly with appropriate chemical potentials. For ex- 
ample, in case of bremsstrahlung, the chemical potential 
for Q and L lines are zero where as chemical potentials 
for R + L, R and P lines are \i q . Similarly, the chem- 
ical potentials for K and K + L lines are \x q for quark 
loop and [i g for gluon loop respectively (see appendix A 
for detail) . The baryo-chemical potential is zero since we 
restrict only to the case of a baryon free plasma. Due 
to the validity of the above factorisation, similar to the 
case of an equilibrated plasma |2^,|2l| , we can show the 
equivalence between two approaches based on field theory 
and kinetic theory. Similar arguments also apply for the 
case of annihilation with scattering due to the following 
identity (given only for quark quark scattering) 

n B {lo)[n B (Zo) + n F (r + Z )] 

x \ n F ( r o) ~ n F (-po)] [n F (k + Zo) - n F (k )} 

= n F (r + l )n F (p )n B (k Q )[l -n F (k + l )}. (34) 

The above identity is quite similar to Eq.(^) except 
that the chemical potential associated with the pre-factor 
n B {qo) is 2[i q . 

Recently, Mustafa et al. |35|] have calculated the pho- 
ton production from a non-equilibrated plasma based on 
the above kinetic theory approach. Assuming a factor- 
ized form of parton distribution functions i.e. n,; = \in e i q 
where n e i q is the distribution function at equilibrium and 
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combining the contributions from quark and gluon scat- 
tering with appropriate spin, colour and flavour statis- 
tics, the rate for the bremsstrahlung production can be 
written as [see Eq.(^2|)]. 

-,X 2 )xlZ 



dR _ 



.1 • /, ,A " Av ' ,1 • B' 



(35) 



where TZ is the equilibrium contribution to photon pro- 
duction, TZ = (A + B)I with J given by Eq.@ to be 
evaluated under classical approximation. The degener- 
acy factors A and B are given by 

4 32 
A=-x2x2 f x2 s x3 c = — ; 
9 ; 3 

5 = 2 S x 8 C = 16. 

The factor 4/9 in the expression for A appears due to the 
assumption l-M^, = 4/9\M\ 2 ^ g where |M|J,_, 9 and 
l-A^lg^g ar e the square of the matrix element for quark- 
quark scattering and quark-gluon scattering respectively 
p7| . Recall that it is possible to combine the quark and 
gluon contribution in a form given by Eq.(|35|) only if 
n F w n B and quantum statistics are ignored in Eq.(^) 
which is true under classical approximation. In this con- 
text, the factorisation in Eq.(p5[) is only approximate. 
This has been the basis of the result used in Ref. [B5j , al- 
though they subsequently use correct equilibrated value 
for TZ obtained from the imaginary part of the photon 
self energy. Note that in Ref. [f35|, instead of 2/5 and 3/5 
as in Eq. (|35|) , these factors are found to be 3/7 and 4/7 
respectively. Apart from this minor discrepancy, the use 
of equilibrium value for TZ is incorrect for the following 
reason even though Eq.(|35|) is approximately correct as 
mentioned before. The factor 1Z contains the square of 
the matrix elements involving the product of the terms 
(R 2 - M^) and (P + L) 2 - in the denominator. Due 
to the presence of two very near by poles, the matrix 
element will be enhancement by an factor of ~ p 2 /M^ 
while the enhancement would have been logarithmic if 
only one of the terms appear in the denominator. There- 
fore, there is an additional enhancement ~ (2 X g + X q )/3 
if the plasma is unsaturated. This additional enhance- 
ment will be compensated to a large extent by the sup- 
pression factor given in Eq. ( |35| ) above particularly when 
the plasma is gluon dominated. In addition to taking 
out the fugacity factors from the distribution functions 
appearing in Eq. ([[2]) , the square of the amplitude should 
also be calculated properly for a non-equilibrated plasma. 
Although a detailed calculation needs to be carried out, 
naively, 1Z should differ from the equilibrium value by 
a factor of ~ A,: 1 when A, « A 9 « 1. A comparison 



with Eq.(35) suggests that the photon production rate 
for bremsstrahlung should have strong dependence on X q 
which (within above approximations) is consistent with 
the results obtained in section III. Based on the similar 
arguments, it can be shown that the AWS photon pro- 
duction will depend only on X 2 . 



D. Suppression and Enhancement mechanism 

Let us now summarize why the suppression and par- 
tial compensation occur when the plasma is unsatu- 
rated. First, it may be easier to understand from the 
kinetic theory arguments. For example, consider the 
case of bremsstrahlung due to quark-quark scattering [see 
Eq.(^||) and Eq.(|33|)]. The suppression is due to the un- 
saturated distribution functions n F (ro + Iq) and n F (ko) 
that appear in the initial states. Similar suppression oc- 
curs in case of quark-gluon scattering except n F (ko) has 
to be replaced by n B (ko). The above suppression is par- 
tially compensated by an additional enhancement factor 
that arise due to the mass effect in the matrix element. 
It may be noted here that the gluon distribution that 
appears in the initial and final states comes through cut- 
ting the effective gluon propagator. In the field theoret- 
ical description, n B (ko) does not appear explicitly, but 
contained in the effective gluon propagator through the 
thermal gluon mass m 2 . Further, notice that the soft 
gluon associated with n B (Iq) and responsible for scatter- 
ing has no role in the above suppression. 

The same thing happens in the field theoretical de- 
scription as well. Since R + L and R correspond to the 
same quark in the initial and final states of scattering, 
we assign same fugacity to both R + L and R lines. This 
argument is also consistent with the kinetic theory for 
the fact that any of the physical processes do not change 
the quark contents in the final states. Therefore, it is 
reasonable to assume the same chemical potential both 
for the initial and final quark lines i.e. to both R + L 
and R lines. From the conservation of potential, it fol- 
lows that L line should have zero chemical potential or 
unit fugacity. Therefore, n B (Iq) still follows the Bose dis- 
tribution. Then, how do we understand the suppression 
due to gluon fugacity which gets compensated by the ad- 
ditional linear enhancement ? Let us consider Eq.feq) 
again. In case of equilibrated plasma, it is shown EG] 
that Eq.(p8|) depends only on the ratio m^/M^. The 
integrals Jt and Jl diverge if the thermal quark mass 
Mqo is switched off. Therefore, the above singularity is 
regularized by this thermal quark mass. In case of non- 
equilibrated plasma, decreases giving enhancement 
in Jt and Jl- However, this enhancement is compen- 
sated by a simultaneous decrease in m 2 . To be more 
explicit, differs from the equilibrium value (g 2 T 2 /6) 
by a factor of (2A g + X q )/3. Similarly, m 2 differs from the 
equilibrium value (4g 2 T 2 /9) by a factor of (3A g + A g )/4 
when Nf = 2. For A g = A g = 1, the ratio m 2 /M^ is 1.33. 



X q <C A g and A s 



« A, 



?! 



Under two extreme conditions 
the above ratios are 1.5 and 1.0 respectively. The cor- 
responding changes in the Jt and Jl values are within 
5% to 10% of the equilibrated values and can be ignored. 
Therefore, it is fair enough to say that the enhancement 
due to the decrease in gets (nearly) compensated by 
the corresponding decrease in m 2 . In other words, for 
a gluon dominated plasma, the collinear enhancement 
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and the suppression due to gluon fugacity cancels out 
leaving Jt and Jl practically unchanged. The only fac- 
tor that suppresses the yield is the p integral in case of 
bremsstrahlung that involves only the quark and anti- 
quark fugacities and the prefactor n B (E) in case of an- 
nihilation with scattering. 



E. One loop results for comparison 

For the completeness and also for comparison, in the 
following, we briefly mention our previous results for ef- 
fective one loop level for the case of a non-equilibrated 
plasma at zero baryon density. The physical processes of 
photon production (annihilation qq — > 97 and Compton 
processes q(q)g — > 5(9)7) in the lowest order of pertur- 
bative expansion (0(aa s )) are obtained by cutting the 
loop diagram as shown in Figure 1. Since the self energy 
is IR divergent in the soft momentum limit, a cut-off pa- 
rameter k 2 is introduced to separate the soft from the 
hard momenta of the intermediate quark. The soft part 
is obtained from a resummed quark propagator and it 
involves a thermal quark mass which acts as an infrared 
cutoff. The hard part is obtained from the relativistic 
kinetic theory from the expression 



2E 



dR 



hard 



N f d 3 p\ d 3 p 2 d 3 p 3 



d 3 q (2tt) 8 
x(l±n 3 (E 3 ))S{p1 



1E\ IE2 2i?3 
P2 - Pa 



" "- <t)Y,\ M 



rii (.Ei) n 2 (E 2 ) 
(36) 



by carrying out integration above the cut-off |l3|Jl6|Ji8[ . 
In the above, 711,2,3 are the parton distribution functions 
with plus sign for annihilation and the minus sign for the 
two Compton processes. The total rate can be obtained 
by adding soft and hard contributions together. The cut- 
off dependence cancels out in the summation. It is also 
found that the total photon rate can be obtained from 
the hard part alone by using the lower limit of integra- 
tion k 2 equal to 2m 2 where the thermal quark mass m 2 
is given in Eq. (gg). Therefore, Eq.(|36"|) can be used to 
estimate the rate of photon production using appropriate 
distribution functions and thermal quark mass for non- 
equilibrium plasma. Although the Juttner functions for 
parton distributions can be used, we restrict to the fac- 
torized form for convenience. However, our conclusions 
are independent of the above choice. Using factorized 
distributions and the identity, 



ni n 2 (l ±n 3 ) 



X 1 X 2 X 3 n e 1 q 
+AiA 2 (l - 



n 2 q {l±nf) 



A 3 )< 



(37) 



the above equation can be broken into two parts 

For the first part, one can use the analytic form that can 

be obtained using Boltzmann distribution for n\ and n 2 

El 



2E 



dR 
d?q~ 



5a a s 
"9^2" 



rj-iZ 



-E/T 



X * X J 2 _ ln (^) _ ] ]:] 



k 2 J 



annihilation 





f\ ET\ 








+ 0.0075 j 







(38) 



C ornpton 

Following Ref. (Taj, the second part can be written as 



2E—\ = 10a as T 2 e~ E/T 



d 3 qj 2 9 7T 
Afl(l-A S K -2-2/3 + 2 In 



4 ET" 
~fcT. 



annihilation 



+ \\ g {\-\)\ 1-2/9 + 2 In 



4 ET\ 



(39) 



C ornpton 



where a, a s are electromagnetic and strong coupling con- 
stants respectively and (3 = 0.577 is the Euler constant. 
In the above, the first term is the contribution from the 
annihilation whereas the second and third terms are due 
to Compton like processes. The total rate is estimated 
by adding Eq.(^8^ and Eq. (|39|) . It may be pointed out 
here that Baier et al. [jl7| have also estimated the rate for 
a non-equilibrated plasma. Although the order of magni- 
tude is same, the above expressions are different from the 
result given in ||l7f due to Boltzmann approximations for 
the initial states. We have also compared the above re- 
sults with the exact numerical calculations using Eq. ( |3^ ) 
and find good agreement. Therefore, we prefer to retain 
the above form for consistency with our previous work 
& 



IV. RESULT AND DISCUSSION 

In the following, we estimate numerically photon pro- 
duction rates for various processes. The crucial aspect of 
the calculation is the numerical integration of the Jt and 
Jl functions which depend sensitively only on the ratio 
m 2 JMl. 
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FIG. 4. The variation of J T and J L with ml/M^ 

Figure || shows the plot of Jt and Jl over a range of mass 
ratios which is of interest in the present study. Although 
Jt and Jl diverge logarithmically in the limit — > 0, 
it has approximately a linear dependence in the above 
region. In case of equilibrated plasma, the above mass 
ratio is about 1.33 (for N — 3 and Nf — 2) and the 
corresponding Jt and Jl values are found to be 1.108 and 
-1.064 respectively |p4| . Note that these values are less 
exactly by a factor of 4 than the values originally reported 
in p"o| , pl"| and used subsequently by many others. Since 
the variation of nig/M^, with fugacity is not significant, 
we also use the above values for Jt and Jl both for 
equilibrated and non-equilibrated plasma. 

Figure |^(a) shows the comparison between one and 
two loop contributions to photon self energy evaluated 
with the corrected values of Jt and Jl- As in J2l]] 
bremsstrahlung dominates in the low momentum region 
whereas AWS dominates in the higher momentum scale. 
Figure ||b shows the photon production rate at a fixed 
temperature T = 0.57 GeV for a chemically equilibrated 
plasma (A s = X q = 1.0). The two loop contribution 
(bremsstrahlung + AWS) competes or even dominates 
over one loop photon production over a wide energy 
range. Further, it is noticed that the bremsstrahlung 
process has strong contribution to photon production be- 
low E ~ lGeV and falls at a faster rate as compared to 
the one loop contribution particularly at higher energy. 
Since the Jt and Jl factors are same, the bremsstrahlung 
and AWS photon productions differ only due to differ- 
ent phase space factors. The p-integral in Eq.(|30|) in- 
volves the quark distribution functions n F (p) — n F (p + E) 
whereas the p-integral in Eq.(^l]) is nearly independent 
of distribution functions for < p < E and has insignifi- 
cant contribution for p > E. Therefore, the phase space 
suppression is stronger for bremsstrahlung as compared 
to the AWS process. The most significant contribution 
to photon production comes from the AWS process at 
higher photon energies. Our results obtained with the 
correct Jt and Jl values are qualitatively in agreement 
with the conclusions drawn from the earlier studies of 
Aurenche et al. 



Next, we consider a chemically unsaturated plasma 
with two different initial conditions at RHIC bombard- 
ing energy. Figure ^(a) corresponds to the initial con- 
ditions T = 0.57 GeV, X g = 0.09 and V = 0.02 as ob- 
tained from HIJING model calculation pa] whereas Fig- 
ure ^(b) is plotted with the initial conditions T = 0.67 
GeV, X g = 0.34 and A ? = 0.064 corresponding to a typi- 
cal SSPC model ||f|. Note that in both cases the initial 
plasma is gluon dominated. A general observation is that 
the AWS contribution is less than the one loop contribu- 
tions and the bremsstrahlung seems to be the dominant 
mechanism of photon production over a wide range of 
energy particularly when the plasma is strongly unsatu- 
rated. The above results can be understood as follows. 
The dominant contributions to photon production at the 
one loop level comes from the term which is linear in 
XgXg [see Eq.© and Eq.©]. Also see Eq.(45) of Ref. 
p7[ . Since the Jt and Jl are insensitive to fugacities, the 
contributions at the two loop level (say) bremsstrahlung 
is suppressed by a factor of X q . 
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FIG. 5. (a) The photon self energy for one loop and two 
loop processes as a function of photon energy, (b) The rate of 
photon production for an equilibrated plasma (A 9 = X q = 1.0) 
at constant temperature T = 0.57 GeV 
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Similarly, the AWS process is suppressed by a factor 
of Xg that arises due to the prefactor although p-integral 
is nearly independent of any distribution functions. As 
a consequence, the AWS process is suppressed strongly 
as compared to the both one loop and bremsstrahlung 
processes. Interestingly, it is the bremsstrahlung which 
dominates the photon production at all energies. This 
is to be contrasted with the equilibrium situation where 
the AWS is the dominant mechanism of photon produc- 
tion at higher energies. Naively, from the kinetic the- 
ory arguments, it is expected that for an unsaturated 
plasma (gluon dominated), the bremsstrahlung and the 
AWS processes will be reduced by a factor of \ q X g and 
Xq\ g respectively. However, the suppression due to X g 
gets compensated by the collinear enhancement which 
goes up by a factor of A~ 1 . These findings are also quite 
intuitive with processes initiated by one quark in the ini- 
tial state being less suppressed than those initiated by 
two quarks. Therefore the suppression factors for the 
AWS, one- loop and the bremsstrahlung processes are A^, 
X q Xg and A^ respectively. 
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FIG. 6. The rate of photon production at (a) T = 0.57 
GeV, A 9 = 0.09, X q = 0.02 corresponds to HIJING model 
and at (b) T = 0.67 GeV, A 9 = 0.34, A g = 0.064 corresponds 
to SSPC model 



V. CONCLUSION 

The effect of chemical potential on photon production 
from a quark gluon plasma has been studied. Since, the 
non-vanishing chemical potential characterizes an unsat- 
urated phase space, the contributions to photon produc- 
tions both at effective one and two loop levels are sup- 
pressed as compared to their equilibrium counterparts. 
The contributions at the effective two loop level i.e. the 
bremsstrahlung and annihilation with scattering (AWS) 
processes are suppressed by a factor of X q and X q respec- 
tively. Interestingly, the above suppressions are found 
to be independent of the gluon fugacity X g . The reduc- 
tion in the photon production rate due to unsaturated 
gluon distribution gets compensated to the large extent 
by the collinear enhancement. This aspect can be under- 
stood from the kinetic theory formalism. However, there 
is no such enhancement for annihilation and Compton 
processes at the one loop level and the rate of photon 
production is suppressed by a factor of ~ A 9 A g . There- 
fore, in case of an unsaturated plasma, the AWS process 
is more suppressed as compared to the one loop contri- 
butions which itself is less by a factor of X g compared 
to the bremsstrahlung process. This is in contrast to 
the equilibrium scenario where AWS dominates the pho- 
ton production followed by one loop and bremsstrahlung 
contributions. In either case, whether the plasma is sat- 
urated or unsaturated, the two loop contributions seem 
to dominate over the one loop processes particularly at 
higher photon energies. 

Further, we would like to mention here that an inher- 
ent assumption which has gone into the above formalism 
is the infinite lifetime of the plasma. As a consequence, 
the photon production rate is independent of time and 
depends only on the photon energy and the temperature 
of the plasma. The consideration of the finite life time 
of the plasma will lead to the time dependent production 
rate which may enhance the photon production further. 
The finite life time effect has been studied in |4j| where 
the plasma is assumed to be both in thermal and chem- 
ical equilibrium. Although the emission rate is a non- 
equilibrium phenomena, the present study is quite differ- 
ent in the sense that the non-equilibrium here refers to 
a chemically unsaturated plasma that evolves with time. 
The basic production rate is still static, but the time de- 
pendence arises due to the hydrodynamical evolution of 
the plasma. Therefore, a meaningful quantity that can 
be compared with the experimental results is the space 
time integrated photon yields from a plasma undergoing 
chemical equilibrium. Such a study is being carried out 
and will be published elsewhere. 
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APPENDIX: RA FORMALISM AT FINITE 
CHEMICAL POTENTIAL 

We generalize the RA formalism [ p6p^ ] appropriate 
for an unsaturated quark gluon plasma (QGP) at finite 
baryon density. Since the QGP is in a thermalized state, 
the parton distributions can be described by the Juttncr 
functions with non-vanishing chemical potential \i. This 
\x can be decomposed as a sum of two components pf 
and fi b where /x c characterizes the unsaturated properties 
and fi b is associated with the finite baryon density of the 
plasma. 

The propagators in RA formalism are 2x2 diagonal 
matrices, constructed from the retarded and advanced 
propagators of the T=0 theory while all the temperature 
dependence appears in the vertices. The propagators 
for fermions and (gauge) bosons, defined on the contour 
characterized by cr, can be written as, 

S F (P) = (tf + M) U [ri] {P) D{P) V [ri] (P) (Al) 
G% V (P) = -g^ U [r >\P) D{P) V [n] (P) (A2) 



where [rj] = B(F) for bosonic (fermionic) 
and D{P) is the diagonal matrix 



propagators 



D(P) = ^ «^ A ; (p) ) (A3) 
with the retarded and advanced propagators given by 



Ar,a(P) 



P 2 - M 2 ± iepo 
The matrices U and V are defined as 

b- 1 



(A4) 



C/M(P) =e f3( - p °- fl) n [r i ] ( Po ) 



c" 1 



(A5) 



(A6) 



where b and c are arbitrary scalar functions of P, r\ = 
±1 for a boson(fermion) and n^po) is the usual Bose- 
Einstcin or Fermi-Dirac distribution defined as 



,['/] 



(Po) 



1 



x = P(po - m) 



(A7) 



(3 being the inverse of temperature and \x is the chemical 
potential as defined above. In the RA formalism, U is 
associated to an outgoing line while V is associated to 
an incoming line. All the temperature dependence which 
is contained in U and V will then appear in the vertices. 



The different types of vertices are calculated depending 
on the momentum flow. Let us consider a vertex with all 



lines incoming as shown in Figure Al 





(a) (b) 

FIG. Al. (a) EM and (b)Strong vertices 

The new vertex function has the form 

-i laf}p (P,Q,R) = -ig M Vf a (P) Vf b (Q) V p F d (R) (A8) 

The Greek indices take the value R or A and the Latin 
indices refer to the 1 (particle) 2 (ghost) of the usual 
formulation of the real time formalism (RTF) so that 
gm = C, <?222 = — Ci where C,—e or g depending on the 
electromagnetic or strong vertex and all other couplings 
being zero. From the definitions above it can be shown 
that, 

l af 3 P (P, 0, R) = ([b(P)} s ^ [b(Q)] 5 e* [b(R)] 5 ** 



x [-c(P)] 5 « A [-c(Q)] hA l-c(R)] SpA 



l-(-l) 



SaR + $pR „-0 L'(, 



(A9) 



with L = p S aR + q 8 f m + r Q 8 pR and L' = [p - 
e(P)pp]S aR + [q -t(Q)pQ}$i3Ft. + [ro -e(R)ii R }S pR . In the 
above, with each momentum P, Q and R, we have intro- 
duced an associated chemical potential e(P)fip, e{Q)^iQ 
and e(R)fi R . The sign function has been introduced to 
ensure that when the momentum reverses, the associated 
chemical potential also changes its sign. This aspect is 
also consistent with the definition of parton distribution 
functions as given in section II. The causality require- 
ment that three particles propagating forward in time 
(or backward in time) can not annihilate into (or be cre- 
ated from) the vacuum demands that jaaa and j RRR 
should vanish. It is immediately clear from Eq. ( |A9| ) that 
J AAA always vanishes. However, the vanishing of ^ RRR 
requires energy and chemical potential conservation. In 
case of finite baryon density, both p c and p, b needs to 
be conserved separately. Therefore, the following set of 
conservation equations are satisfied when "f RRR = 



Po - 
p c P 



q +r = 0; 
H + Mij = 0; 



fjp + f i b R = 



(A10) 



Note that in the above the baryo-chemical potential for 
photon or gluon has been set to zero. Next, we c onsider 
a crossing fermion line as shown in Figures A2(a) and 
|A2|(b) with the conservation laws, 
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P + Q = R; p c P + pf Q = p R ; p b P = p b R (All) 
The vertex function can be evaluated from 

-i^ P {P,Q;R) = -ig abd V^ a (P) vf b (Q) u^(R) 

(A12) 

Using the definition of U and V, the abov e equation can 
be written similar way as that of Eq. ( |A9[ ) given by, 

7a/3;p (P,Q;R) = C[b(P)} s " R [b(Q)} 5 ^[b(R)}- 5 » R 



x[-c(P)e 



-<rp -\5 aA f 



n F (r )[e^ ro -^] s " B (-1)'"* + (-1) 5 " R (A13) 

where, Pq = [p - e(P)fi P ]S aR + [q - e(Q)p Q }S l3R - [r - 
e(R)fj, R ]5 p A- 




FIG. A2. (a) and (b) are the QED and QCD vertices 
for crossing fermion line (c) and (d) are the QED and QCD 
vertices for crossing boson line 



By comparing with Eq.(A9), we can derive the rela- 
tions 



7a/3;R(P, Qi R) 



e~ gr ° n F (-r ) 
b(R)c(-R) 



1 af )A{P,Q,-R) 



(A14) 



'Y a f3;A(P,Q;R) = 

The choice 



n F (r ) 



b(-R)c{R) 



l a0R (P,Q,-R) (A15) 



b(-R)c(R) = -n F (r Q )e° r ° (A16) 

gives the crossing relation for fermion 

7ct/3;p (F, Q; R) = laf3 p(P, Q, -R) (A17) 

where p — A, R is the conjugate index of p — R, A. The 
crossing property of boson [see Figures |A2|(c) and A2 (d)] 



can also be derived in a similar way except the conserva- 
tion 

P + R = Q; pp + pf R = p c Q ; p b P + p b R = (A18) 

sho uld b e followed. The replacement of — np(ro) in 
Eq.(Al6) with n B (qo) leads to the equation for boson 



b(-Q)c(Q)=n B (q )e°i°. 



(A19) 



It needs to be stressed here that the factor n B (go) rep- 
resents the boson distribution function of energy q and 
with appropriate chemical potential pq. This distribu- 
tion is not to be associated with the emitted photon of 
energy qo which has no c hem ical potential and distribu- 
tion fu nctio n. Using Eq.([A9| ) with the conditions given 
by Eq.( Al6| ) and Eq.( A19) and also with the choice 6=1, 
all the required vertices can be calculated. For example, 
we consider Figure A3 which contributes to the physical 
process annihilation with scattering. 




FIG. A3. Self energy which contribute to AWS process 



The Eq.(A9) can be simplified for e using appro- 



priate conservation laws to get 



^AAB. 



(R,P,-Q) 



e n F (r )n F (p ) 



(A20) 



which has same form what one would have expected 
for an equilibrated case except the distribution func- 
tions which now contain appropriate chemical poten- 
tial. Assuming pp — p c R — p c q (quark chemical po- 
tential) and \p P \ = \p b R \ = p b , the conservation laws, 
Q = P + R; Pq = pp + p c R and p b P + p b R = suggest 
that the boson line should have a total chemical potential 
2p c q while the total chemical potential for P and R lines 
are p q + p b and pf — p b . (Note that in this topology, 
the chemical potential p R and pp are associated with 
a quark and anti-quark respectively. Both should have 
same chemical potential p c and opposite baryo-chemical 
potential p b ). Using the distri butio n function with the 
above chemical potentials, Eq.(|A20|) can also be written 



„AAR 



(R, P, -Q) = e[n F (r Q ) - n F (-p )} (A21) 



Similarly, we can write other vertices 

g AAR (R + L, —L, -R) = g[n B (l ) + n F (r + l )] 
e^R{Q,-P + L,-R-L) 



ARR 



(Q, -P + L, -R - L) = -e 



„ARR, 



(P - L, -P, L) = q {P - L, -P, L) = -g 

(A22) 

Since the chemical potential associated with R+L and R 
are same (p = p q +p b ), the chemical potential for L line is 
zero. The above diagram contributes to bremsstrahlung 
when P reverses it's direction. In this topology, the total 
chemical potential for R, P and Q lines are p q + p b , p q + 
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H h q and respectively. The corresponding four vertices 
are 

e AAR (R, -P, -Q) = e[n F (r ) - n F (p )] 
g AAR (R + L, —L, -R) = g[n B (l ) + n F (r + l )\ 

e ARR (Q, P + L, —R — L) 



e ARR (Q^ P + L, -R - L) 



= —e 

,.ARR 



(-P-L,P,L) 



ARR 



(-P-L,P,L) 



-g 

(A23) 
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